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ABSTRACT 

The apparent shapes of spiral galaxies in the 2-Micron All Sky Survey Large Galaxy 
Atlas are used to constrain the intrinsic shapes of their disks. When the distribution 
of apparent axis ratios is estimated using a nonparametric kernel method, the shape 
distribution is inconsistent with axisymmetry at the 90% confidence level in the B band 
and at the 99% confidence level in the Ks band. If spirals are subdivided by Hubble 
type, the late-type spirals (Sc and later) are consistent with axisymmetry, while the 
earlier spirals are strongly inconsistent with axisymmetry. The distribution of disk 
ellipticity can be fitted adequately with either a Gaussian or a lognormal distribution. 
The best fits for the late spirals imply a median ellipticity of e ~ 0.07 in the B band 
and e ~ 0.02 in the Kg band. For the earlier spirals, the best fits imply a median 
ellipticity of e « 0.18 in the B band and e « 0.30 in the Kg band. The observed scatter 
in the Tully-Fisher relation, for both late and early spirals, is consistent with the disk 
ellipticity measured in the B band. This indicates that excluding spirals of Hubble type 
earlier than Sc will minimize the intrinsic scatter in the Tully-Fisher relation used as a 
distance indicator. 

Subject headings: galaxies: spiral - galaxies; photometry - galaxies: fundamental 
parameters 


1. Introduction 

The intrinsic shapes of disks in spiral galaxies, although difficult to determine, are useful to 
know. In particular, the ellipticity of a disk in its outer region reflects the nonaxisymmetry of 
the potential in which it exists, and thus provides information about the shape of dark matter 
halos. Non-circular disks produce scatter in the Tully-Fisher relation between rotation velocity 
and absolute magnitude (Tully & Fisher 1977; Franx & de Zeeuw 1992); therefore, knowledge of 
the disk ellipticity helps in understanding the origin of scatter in the Tully-Fisher relation. In 
particular, the usefulness of the Tully-Fisher relation as a distance indicator would be enhanced for 
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a population of disk galaxies with small intrinsic ellipticity. If the shape of a disk is approximated 
as an ellipsoid with principal axes of length a > b > c, the disk shape can be described by two 
parameters: the dimensionless disk thickness, 'y = cja, and the disk ellipticity, e = 1 — 6 /a. The 
simplest approximation, that disks are infinitesimally thin and perfectly circular (7 = 0, e = 0), is 
occasionally useful; for many purposes, however, it is useful to have a more accurate estimate of 
the shapes of disks. 

The thickness of disks can be determined from catalogs of edge-on spiral galaxies, such as 
the Flat Galaxy Catalogue (Karachentsev et al. 1993) and the Revised Flat Galaxy Catalogue 
(Karachentsev et al. 1999). Some superthin galaxies have axis ratios as small as 7 ~ 0.05 (Goad & 
Roberts 1981). The study of Guthrie (1992), using images of edge-on spirals from blue Palomar Sky 
Survey plates, found that \ogq had a mean value of —0.95, where q is the apparent axis ratio. The 
mean disk thickness depends on the spirals’ morphological type (Bottinelli et al. 1983); Guthrie 
(1992) found that the mean value of log q ranges from —0.7 for edge-on Sa galaxies to — 1.1 for edge- 
on Sd galaxies. The disk thickness of an individual spiral galaxy also depends on the wavelength at 
which it is observed. Edge-on spirals have strong vertical color gradients, becoming much redder 
with increasing height above the midplane (Dalcanton &: Bernstein 2002). This gradient reflects 
the transition from the thin disk, with its young blue stellar population, to the thick disk, with 
its older redder population. As a consequence, q of an edge-on galaxy becomes larger at longer 
wavelengths. For example, the galaxies in the Dalcanton & Bernstein (2000) sample were chosen 
to have q < 0.125 at blue wavelengths; in deep R-band images, the typical axis ratio increases 
to O' ~ 0.23 (Dalcanton &: Bernstein 2002). Similarly, Mitronova et al. (2003) found that edge-on 
spirals that have g ~ 0.17 in blue images have q ~ 0.35 in two-micron images. 

Estimates of the distribution of disk ellipticity e can be made from the distribution of apparent 
axis ratios q for a large population of randomly oriented spiral galaxies. The signature of non¬ 
zero ellipticity, in this case, is a scarcity of nearly circular (q ~ 1) spiral galaxies. Binney &: de 
Vaucouleurs (1981), using spiral galaxies from the Second Reference Gatalogue of Bright Galaxies 
(de Vaucouleurs, de Vaucouleurs, & Gorwin 1976), concluded that late-type spirals were better 
fitted by slightly elliptical disks (e ~ 0.1) than by perfectly circular disks. More recent studies 
(Grosbpl 1985; Lambas, Maddox, & Loveday 1992; Fasano et al. 1993; Alam &: Ryden 2002; Ryden 
2004) have confirmed the ellipticity of disks, and have permitted estimates of the distribution 
function for intrinsic ellipticity. The distribution of e is commonly fitted with a half-Gaussian 
peaking at e = 0. Values of erg, in this case, range from cje = 0.12 (Lambas, Maddox, & Loveday 
1992), through = 0.13 (Fasano et al. 1993), to = 0.21 (Ryden 2004). 

Kinematic studies can also provide an estimate of the ellipticity of the gravitational potential 
in the disk plane. In an elliptical potential with small ellipticity e#, the closed stellar orbits will 
themselves be ellipses (Binney 1978). When the resulting elliptical stellar disk is seen in projection, 
the isophotal principal axes will be misaligned with the kinematic principal axes. Because of the 
misalignment, there will be a velocity gradient along the isophotal minor axis proportional to 
sin 2(/>, where cj) is the azimuthal viewing angle relative to the long axis of the potential (Franx 
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&: de Zeeuw 1992). Measurement of the kinematic misalignment for a sample of 9 disk galaxies 
yielded an average ellipticity of e<j> ~ 0.044 (Schoenmakers, Franx, & de Zeeuw 1997; Schoenmakers 
1998). A similar analysis of gas kinematics in 7 disk galaxies by Beauvais &: Bothun (1999) found 
an average ellipticity of e ~ 0.059 for the gas orbits. Andersen &: Bershady (2002), using the 
method of Andersen et al. (2001), combined two-dimensional kinematic and photometric data for 
a sample of 28 nearly face-on disk galaxies; they found an average halo ellipticity of e$ « 0.054 and 
an average disk ellipticity of e ~ 0.076 for their sample. Taking into account the selection criteria 
of Andersen &: Bershady (2002), the implied median disk ellipticity is e ~ 0.10 (Ryden 2004). 

Much of what we know about the ellipticity of disks thus comes from looking at galaxies with 
low inclination; much of what we know about disk thickness comes from looking at galaxies with high 
inclination. In principle, both ellipticity and thickness can be determined by looking at a sample 
of disk galaxies with random inclinations (Fasano et al. 1993). However, at visible wavelengths, 
thanks to the effects of dust, it’s difficult to create a sample in which the disk inclinations are 
completely random. Because disk galaxies are not fully transparent, a flux-limited sample will be 
deficient in nearly edge-on galaxies (Huizinga & van Albada 1992). Because disk galaxies are not 
fully opaque, an angular diameter-limited sample will have an excess of nearly edge-on galaxies 
(Huizinga & van Albada 1992; Masters, Giovanelli, & Haynes 2003). 

To minimize the effects of dust, I will be examining the apparent shapes of disk galaxies in the 
Two Micron All-Sky Survey (2MASS). In §2, I extract a randomly inclined sample of spiral galaxies 
from the 2MASS Large Galaxy Atlas. In §3, I make both nonparametric and parametric estimates 
of /(y), the distribution of disk thickness, assuming the disks are axisymmetric. Separate estimates 
are made for early-type spirals and for late-type spirals. In §4, I make nonparametric estimates of 
/(e), the distribution of disk ellipticity for nonaxisymmetric disks. Both Gaussian and lognormal 
functions are found to give satisfactory fits to /(e). In §5, I examine the influence of the deduced 
disk ellipticity on the scatter of the Tully-Fisher relation, and consider the possible origins of disk 
ellipticity for early-type and late-type spirals. 


2. The Galaxy Sample 

The 2MASS Large Galaxy Atlas (LGA) is a catalog of galaxies with large angular size as seen 
in the near-infrared (Jarrett et al. 2003). The LGA provides an essentially complete sample of 
galaxies with 7220 ^ 60'^, where 7220 is the semimajor axis of the Kg band 20 mag arcsec”^ isophote. 
Observations in the Kg band, with an effective wavelength of 2.2/im, have two advantages for 
the study of galaxy morphology. First, the effects of dust are minimized. Second, the light in 
the Kg band emphasizes the smoother distribution of population H stars rather than the patchier 
distribution of population I stars, which tend to lie along spiral arms. The LGA galaxy images 
consist of individual 2MASS co-add images mosaicked together. The typical la background noise 
in a single Kg co-add image is ~ 20mag arcsec”^; the point spread function FWHM varies from 
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2 to 3 arcseconds, depending on seeing (Jarrett et al. 2003). The on-line Large Galaxy Atlas^, 
version 2.0, contains 620 galaxies. To create a sample of spiral galaxies, I selected those galaxies 
with morphological classification SA, SAB, SB, or S. Galaxies labeled as ‘pec’ were excluded, as 
were SO galaxies. The resulting sample contains 383 spiral galaxies. 

The apparent shape of each galaxy in the Kg band is estimated by Jarrett et al. (2003) 
from the axis ratio of the 3 cj intensity isophote, corresponding to a surface brightness of fix ~ 
18.8 mag arcsec“^; I call this axis ratio qk- For purposes of comparison, the apparent axis ratio 
qb in the B band was estimated by taking the axis ratio of the = 25 mag arcsec“^ isophote 
from the LEDA database. Because the 2MASS Large Galaxy Atlas is an angular diameter-limited 
sample and spiral galaxies are largely transparent in the Kg band, the Large Galaxy Atlas contains 
an excess of high-inclination disks with small qx- For instance, of the 30 LGA galaxies with the 
largest value of R 20 , 10 are spiral galaxies with qx < 0.33. 

To account for inclination effects, we must estimate what the isophotal radius of each galaxy 
would be if we saw it face-on. Suppose, to begin, that the isoluminosity surfaces of the galaxy are 
adequately approximated as similar, concentric, coaxial ellipsoids, given by the formula 

m^ = x^ + yV/9' + , (1) 

where 1 > /9 > 7 . The galaxy is viewed from an angle {9, (p), where 9 is the polar angle measured 
from the positive z axis and (p is the azimuthal angle measured from the positive x axis. Assuming 
perfect transparency, the isophotes of the galaxy will be ellipses, given by the formula 

= DX"^ + EY"^ . (2) 

If I define, following Binney (1978), 

A{9, (p) = cos^ 6*(/3^ sin^ cp + cos^ cp) + 7 ^ sin^ 9 , (3) 

B{9,(p) = {I — p'^) cos 9 sm2cp , (4) 

C{(p) = sin^ (p + 0^ cos^ (p , (5) 

and 

f{9, (p) = 7 ^ sin^ 9{0^ cos^ cp + sin^ cp) -|- 0^ cos^ 9 , ( 6 ) 

then 

D{9, cP) = ^^[A + C- V(A-C)2 + B2] (7) 

and 

E{9, cP) = ^\A + C+^{A- 0)2 + B2] ( 8 ) 


^See http://irsa.ipac.caltech.edu/applications/2MASS/LGA/ 
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Jarrett et al. (2003) find that LGA galaxies are well fitted, outside their central nucleus or core, by 
a SctsIc profile. 


/i(X) = Aid exp 


(9) 


The functional form adopted by Jarrett et al. (2003), in which a is a scale radius, is readily converted 
to the more usual form, 


Ai(A) = Ai(i?e) exp [-hn[{X/Ref/^ - 1]) , 


( 10 ) 


in which is the effective radius, or half-light radius. The conversion from scale length a to 
effective radius Re is Re = {bn)^a, where an adequate approximation to bn is (Ciotti &: Bertin 
1999 ) 

14 46 

^n-2n-- + ^ + -^. (11) 


For the spiral galaxies in the LGA, the best fitting profiles are usually close to exponential; 90% 
of the spiral galaxies have 0.7 < n < 1.4 in the Kg band. In this range of n, the effective radius 
ranges from Re = l.OGa (when n = 0.7) to Re = 3.55a (when n = 1.4). 

If a galaxy has a surface brightness profile given by equation (9), then if it were viewed face-on 
(that is, from along the 0 = 0 axis) its profile would be 


Ai(A)=/r2exp f-(A/a°)V^ 


( 12 ) 


Here, a° is the scale length along the apparent long axis and At^ is the central surface brightness, 
as seen from the angle 0 = 0. Let R^q be the value of X at which the surface brightness falls to 
the specified isophotal brightness /j, 2 o = 20 mag arcsec“^. For the profile of equation (12), 

r^2o = ■ (13) 

If we actually view the disk galaxy from an angle {0,(p), the central surface brightness of the disk 
is 

^ei{0, 0) = fi2PVDiO,<P)Ei0,cj2) (14) 

and the scale length along the apparent long axis is 

a{0, (/)) = a°/■\/D{0,(j)) . (15) 

The isophotal radius R 20 , at which the surface brightness on the long axis falls to ^ 20 , is 

R2o{0,(p) = a[ln{fid/0-2o)r • (16) 


The apparent axis ratio of the disk, seen from the angle {0,4>), is g = sjD/E. Thus, the face-on 
isophotal radius, i? 20 > related to the observed values of R 2 q{ 0 , 4 >) and a{0,4>) by the relation 


= ^/dR 


20 


l + (a/i22o)'/"ln(^j 


(17) 
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Although R 20 , a, and q can be measured for any particular galaxy seen in projection, /3 and D{6,4>) 
cannot. However, if the disk is axisymmetric, then (5 = D = 1, the scale length a is independent 
of viewing angle, and the face-on isophotal radius can be deduced from the observed properties: 


-^20 ~ -^20 


1 -I- In O' 


(18) 


For all n > 0, an axisymmetric disk appears smallest when seen face-on — -^ 20 )- 


3. Shape of Axisymmetric Disks 

If we assume that spiral galaxies are axisymmetric, then the analysis of the preceding section 
suggests a way of assembling an inclination-independent subsample of disk galaxies from the 2MASS 
Large Galaxy Atlas. If the LGA contains all galaxies with R 20 > 60'', then, assuming axisymmetry, 
it will contain all galaxies with > ^0". Thus, from the tabulated values for R 20 ) n, and q 
(Jarrett et al. 2003), I compute ii^o each spiral galaxy, using equation (18). If R^q > 60", the 
spiral is included in the inclination-independent snbsample. 

In this way, I find that 193 spiral galaxies in the 2MASS Large Galaxy Atlas have R^q > 60". 
Of these galaxies, 130 have a Hubble type earlier than Sc (14 SO/a, 6 Sa, 28 Sab, 40 Sb, and 42 
Sbc); I’ll call these galaxies the “early spirals”. The remaining 63 spirals have a Hubble type of Sc 
or later (36 Sc, 18 Scd, 8 Sd, and 1 Sm); I’ll call these the “late spirals”. If I had not done the 
inclination correction to find R 20 , I would have found 161 early spirals and 117 late spirals with 
R 20 > 60". 

The apparent axis ratios qx and qB for the 193 LGA spirals in the inclination-corrected sample 
are shown in Figure 1. Spirals with qs < 0.3 are generally less flat in the Kg band than in the 
B band. It can also been seen that extremely flat spirals, with {qs + qK)/"^ < 0.2, and nearly 
round spirals, with {qs + qK)l‘^ > 0.9, are more likely to be late types (squares in Figure 1) than 
early types (triangles in Figure 1). Intermediate axis ratios, however, with (g^ -|- qK)l‘l ~ 0.5, are 
strongly dominated by early spirals. When galaxies with detectable visible-light bars (types SAB 
and SB) are compared to galaxies of type S and SA, no significant difference is found in their axis 
ratio distributions; P > 0.1 as measured both by a Kolmogorov-Smirnov test and by a test, in 
both bands. 

On a cautionary note, it should be pointed out that qx^ the axis ratio of the 3cr isophote, doesn’t 
measure the shape at the same physical distance from each galaxy’s center. Because looking at 
a disk edge-on enhances the surface brightness, the distance Pso-, the semimajor axis of the 3cr 
isophote, will be larger when a galaxy is seen edge-on than when it is seen face-on. Thus, qx is 
measuring the shape at larger radii, on average, for galaxies with small qx than for galaxies with 
large qx- To illustrate the trend, Fignre 2 shows R^a, in units of the effective radius. Re = {bn)^a, 
as a function of Inqx- On average, R^a/Re tends to increase as qx decreases. The best linear fit 
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to the data for late spirals, shown as the dashed line in Figure 2, is 

^ = 0.976 -0.617 In , (19) 

with a correlation coefficient r = —0.64. The best fit to the data for early spirals, shown as the 
dotted line in Figure 2, is 

Rn 

^ = 1.235 -0.670 In g/c , (20) 

Re 

with a correlation coefficient r = —0.56. Although the scatter in the relation is large, the 3cr 
isophote tends to be closer to the galaxy’s center for nearly face-on axisymmetric disks than for 
nearly edge-on disks. If the isoluminosity surfaces of a disk galaxy were perfectly similar oblate 
spheroids, it wouldn’t matter where qx is measured. However, for early-type, big-bulged spirals, 
the value of qk measured for galaxies of low inclination can be significantly affected by light from 
the galaxy’s central bulge. A similar relation holds true between i ?25 and qb, the apparent axis ratio 
of the rb = 25 mag arcsec”^. However, since R 25 in the B band is typically larger than i? 3 o- in the 
Kg band, the value of qb is usually not strongly influenced by bulge light. For the early spirals, the 
mean and standard deviation of R 2 b/R‘ir 7 are 1.90±0.62; for the late spirals, R 2 b/R^a = 2.45±1.04. 
(For two early spirals at low galactic latitude, WKK 4748 and Maffei 2, R 2 ^/Rz (7 ^ 0.5.) 


3.1. Nonparametric fits 


The distributions of apparent axis ratios, /{qk) and /{qb), were estimated for the galaxies 
in the sample using a nonparametric kernel technique (Vio et al. 1994; Tremblay &: Merritt 1995; 
Ryden 1996). Estimating a function / by the use of kernels gives a smooth estimate function / 
with no a priori restrictions on its shape, unlike a parametric fit to /. Given a sample of n axis 
ratios, qi, q 2 , ■ ■ ■, qn, the kernel estimate of the frequency distribution f{q) is 


m = 



( 21 ) 


where K(x) is the kernel function, normalized so that its integral over all x is equal to one. To ensure 
that the estimate / is smoothly differentiable, I adopt a Gaussian kernel, K{x) oc exp(—x^/2). The 
parameter hi in equation (21) is the kernel width, determined using the standard adaptive two- 
stage estimator of Abramson (1982). In the two-stage estimator, the first stage consists of making 
an initial estimate /o using a fixed kernel width h. The initial kernel width is given by the formula 
h = 0.9An“*^'^, with A = min(cr, <54/1.34), where a is the standard deviation of the axis ratios and 
<54 is the interquartile range. For samples that are not extremely skewed, this formula minimizes 
the mean square error of the estimate (Silverman 1986; Vio et al. 1994). The second stage consists 
of making the final estimate /, using at each data point q^ the kernel width 


G 

foiQi) 


hi = h 


1/2 


( 22 ) 








where G is the geometric mean of /o over all values of qi. 

The estimated distribution of axis ratios for the 193 spirals in the inclination-corrected sample 
is shown in the upper panel of Figure 3. The solid red line indicates /{qk), while the solid blue 
line indicates f{qB)- The dotted lines indicate the 80% error intervals estimated by bootstrap 
resampling. In the bootstrap analysis, I randomly selected n = 193 data points, with substitution, 
from the original sample of 193 axis ratios; a new estimate / was then created from the bootstrapped 
data. After doing 5000 bootstrap estimates, I determined the 80% error interval; in Figure 3, 10% 
of the bootstrap estimates lie above the upper dotted line, and 10% lie below the lower dotted line. 
If the spiral galaxies were infinitesimally thin, perfectly circular disks, the distribution function 
would be / = 1 for 0 < g < 1; the deviations of /{qx) and /{qs) from this form reveal that the 
spiral galaxies in our sample are neither infinitesimally thin nor perfectly circular. In the Kg band, 
galaxies are slightly rounder on average than in the B band: (qx) = 0.538 versus {qs) = 0.525. 
In addition, there is a slightly greater spread in shapes in the B band: (T{qx) = 0.212 versus 
a{qB) = 0.234. 

Given the assumption of oblateness, it is straightforward to do a numerical inversion of f{q), the 
estimated distribution of apparent axis ratios, to find N{'y), the estimated distribution of intrinsic 
axis ratios. The relation between f{q) and N{'y) is 

Hq) = [ P{qh)N{j)dj , (23) 

Jo 

where P{q\'y)dq is the probability that a galaxy with a short-to-long axis ratio 7 has an apparent 
axis ratio in the range q —> q+dq when seen from a randomly chosen viewing angle. For a population 
of oblate spheroids (Sandage, Freeman, & Stokes 1970), 

= (1 _ ^2)1/2(^2 _ ^2)1/2 

if 7 < O' < 1, and P{q\'y) = 0 otherwise. Equation (23) is a Volterra equation of the first kind; in 
its discretized form, it can be inverted to find A by a process of forward substitution (see Vincent 
& Ryden (2005) for details). 

The bottom panel of Figure 3 shows the distributions of intrinsic axis ratios in the Kg band 
(red line) and the B band (blue line). The solid line, in each case, is the distribution found by 
inverting equation (23). Because iV is a deconvolution of /, any noise present in / will be amplified 
in N] to produce an acceptably smooth estimate of N, it is thus necessary to use a value of the 
initial kernel width h that is larger than the optimal value for computing / (Tremblay &: Merritt 
1995). For computing N, I used a kernel width h that was 20% greater than the value used to 
compute /. 

The deduced distributions N show that the most probable value of 7 is 'yx = 0.25 in the Kg 
band and 7^ = 0.17 in the B band. A striking aspect of the distributions N seen in Figure 3 is that 
they are negative for large values of 7. In the Kg band, N dips below zero for 'yx > 0.58; in the 
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B band, N is below zero for 7 ^ > 0.67. This is a physical absurdity, implying a negative number 
of nearly spherical galaxies. The 80% confidence band, delineated by the dotted lines in the lower 
panel of Figure 3, falls below zero, in both bands, for 7 > 0.7. This means that fewer than 10% 
of the bootstrap resamplings yield non-negative estimates in this range of 7 . The hypothesis that 
the spiral galaxies are a population of randomly oriented oblate spheroids can thus be rejected at 
the 90% (one-sided) confidence level. If the 98% confidence band is laid out in a similar manner, 
with 1 % of the bootstrap resamplings lying above the band and 1 % lying below, it is found that 
in the Ks band, the confidence band dips below zero for large 7 . Thus, in the Kg band, the oblate 
hypothesis can be rejected at the 99% (one-sided) confidence level. In the B band, where the shape 
is measured farther from the galaxies’ centers, the oblate hypothesis cannot be rejected at the 99% 
confidence level. 

For the complete inclination-adjusted sample, the hypothesis that the galaxies are perfectly 
axisymmetric can be rejected at a fairly strong confidence level, particularly in the near infrared. 
However, Figure 1 reveals that the distribution of apparent shapes is different for late spirals and 
for early spirals. Thus, it is worthwhile to look at f{q) and the derived N{^) for late and early 
spirals separately. The upper panel of Figure 4 shows /{qk) and /{qb) for the n = 63 late spirals in 
the sample. The distribution of apparent axis ratios in the Kg band, indicated by the red lines, is 
similar to the distribution in the B band, indicated by the blue lines. In both cases, the distribution 
/ is nearly constant over a wide range of q. The deduced distributions of intrinsic axis ratios, N{'^k) 
and N('yB), are indicated by the red and blue lines respectively in the bottom panel of Figure 4. 
In both bands, the distribution N never falls below zero; thus, the oblate hypothesis is perfectly 
acceptable, both for the B band data and for the Kg band data. If the oblate hypothesis is accepted, 
the average disk thickness is smaller in the B band: {jb) = 0.12 as compared to {jk) = 0.19. In 
both bands, the width of the function N is comparable to the kernel width h used, indicating that 
the true spread in 75 and niay be less than is shown in the smoothed distributions of Figure 4. 

The apparent shapes of early spirals, shown in the upper panel of Figure 5, are strikingly 
different from the apparent shapes of late spirals. For early spirals, the distribution / is not 
uniform over g, but instead is peaked at g ~ 0.5, with a marked scarcity of galaxies with g > 0.8. 
Such a nonuniform distribution is typical of a population of triaxial ellipsoids. The distribution of 
qx (red lines in Figure 5) is particularly strongly peaked at g ~ 0.5, which may be attributed, in 
part, to a contribution of light from nonaxisymmetric bulges. However, even the distribution of 
gs, measured at a much larger radius within each galaxy, is inconsistent with axisymmetry. The 
deduced distributions of intrinsic axis ratios, shown in the lower panel of Figure 5, demonstrate how 
grossly inconsistent the observed distributions of qx and qx are with the hypothesis of randomly 
oriented oblate spheroids. The oblate hypothesis can, in fact, be ruled out at the 99% confidence 
level for both the Kg band isophotal shapes and the B band isophotal shapes. Thus, the deviations 
from axisymmetry in the full sample of spiral galaxies can be attributed to the early spirals in the 
sample; the late spirals are perfectly consistent with axisymmetry. 
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3.2. Parametric fits 


For comparison with previous studies, it is useful to make parametric estimates of N('y), in 
addition to the nonparametric kernel estimates of the previous subsection. A common convention 
(Huizinga & van Albada 1992; Fasano et al. 1993) is to fit the distribution of intrinsic axis ratios 
with a Gaussian: 


N{'y)d'y oc exp 


(7 - 

2a2 


dj 


(25) 


for 0 < 7 < 1, with = 0 elsewhere. Inspection of the lower panel of Figure 4 indicates that this 
should provide a reasonable fit for the late spirals, at least. 


The best-fitting values of and cr-,, were determined by a fit to the binned distribution 
of qk and qb- In each band, the n values of q were divided into nbin = 12 bins of equal width. 
After choosing particular values for and I randomly chose n values of 7 , drawn from the 
distribution of equation (25). After randomly selecting a viewing angle for each mock disk, I 
computed its apparent axis ratio 5 ( 7 , 9, 4>) from the formulas in §2. The model axis ratios were 
then binned in the same way as the observed axis ratios. After repeating this procedure 8000 times 
for a given pair, I calculated the mean and standard deviation of the expected number of 

model galaxies in each bin, and computed a probability for that particular pair of parameters. 


Figure 6 shows the isoprobability contours in the two-dimensional parameter space, where the 
probability is measured by a fit to the binned data. The fit to the Kg band data is shown by 
the red lines; the fit to the B band data is shown by the blue lines. In each band, the best fitting 
values of and cr-,, are indicated by an asterisk. 

In the top panel of Figure 6 , the fit is made to the complete inclination-adjusted sample of n = 
193 spirals. In the Kg band, the best fit in this region of parameter space is = (0.26,0.09), 

which yields a x^ probability P ^ 0.20. In the B band, the best fit is = (0.17,0.03), 

which yields P ~ 0.17. For these fits, the biggest contributions to x^ comes from the two bins 
with q > 0.83, for which this axisymmetric model predicts more galaxies than are observed. In 
the Kg band, spiral galaxies are not only thicker, on average, than in the B band, but they also 
have a greater spread in thicknesses. The typical 7 I find for spirals in the B band, = 0.17, is 
greater than the value of fx-y = 0.13 found by Fasano et al. (1993) for a volume-limited sample of 
766 spirals, using apparent axis ratios taken from the Third Reference Catalogue of Bright Galaxies 
(de Vaucouleurs et al. 1991). As seen in the upper panel of Figure 6 , the best values of {fx-y,a^) 
found by Fasano et al. (1993) lie just inside the P = 0.01 contour for the 2MASS LGA spirals. 

The middle panel of Figure 6 shows the fit for the late spirals only (Hubble type Sc and later). 
The n = 63 galaxies in this subsample do not place a strong constraint on and the P = 0.5 
contour, indicated by the heavy solid line, contains a large swath of parameter space, particularly 
in the Kg band. The best fit in the Kg band is {ix-y,a^) = (0.18,0.06), which yields a x^ probability 
P Ri 0.91. In the B band, the best fit is (/r.^,cr-y) = (0.14,0.02), yielding P Ri 0.80. Once again, as 
for the complete sample, the Kg band data indicates both a larger average thickness and a greater 
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spread in thicknesses than does the B data. However, in both bands, the data are consistent with 
a^y = 0. An excellent fit results if all the late spirals are assumed to have an intrinsic axis ratio 7 
equal to, or slightly smaller than, the smallest apparent axis ratio q in the subsample {qk = 0.12 for 
NGC 4302, and qb = 0.12 for NGC 5907). The best fit found by Huizinga & van Albada (1992) for 
a sample of Sc galaxies, {ny, ay) = (0.16,0.02), and by Fasano et al. (1993) for a sample of galaxies 
of type 4.5 < t < 7.0, {fiy,ay) = (0.10,0.03), are both statistically consistent with the results for 
late-type LGA 2MASS spirals in the B band; see the middle panel of Figure 6 . 

The bottom panel of Figure 6 shows the results of fitting a Gaussian distribution for 7 to 
the early spirals (Hubble type Sbc and earlier). The best fits are not as good as the excellent 
fits found for the late spirals. The best fit in the Kg band is {ny,ay) = (0.29,0.08), for which 
the probability is P ^ 0.06. The best fit in the B band is {ny,ay) = (0.20,0.04), yielding 
P ~ 0.08. The relatively poor fits are the result of demanding axisymmetry; the Gaussian form of 
equation (25) is not primarily to blame, since the nonparametric kernel estimate revealed that no 
distribution N{'y) > 0 gives a good fit to the early spirals. As with the late spirals, the average 
thickness and the spread of thicknesses is greater in the Kg band than in the B band. The greater 
disk thickness in the infrared is more pronounced for the early spirals than for the late spirals. For 
early spirals, the best-fitting ^y in the Kg band is 50% greater than in the B band; for late spirals, 
the best-fitting fiy in the Kg band is only 20% greater than in the B band. 


4. Shape of Nonaxisymmetric Disks 

Given that a distribution of axisymmetric (e = 0) disks does not give a good fit to early spirals, 
it would be useful to know which distribution of nonaxisymmetric disks gives the best fit to the 
early spirals (and to the late spirals, for that matter). A fundamental difficulty is that even if f{q) 
were known perfectly, it would not uniquely determine N{'y,e). However, the joint distribution of 
the disk thickness 7 and the disk ellipticity e is not beyond all conjecture. I will use the approach 
of assuming a parameterized functional form for the distribution N{'y,€), and finding which values 
of the parameters give the best fit to the observed distributions of axis ratios. This approach yields 
permissible fits to N{'y,e), but these fits, it should be remembered, are not unique solutions. 


4.1. Gaussian parametric fits 


I will assume that N{'y,e) = Ny{'y) x A'e(e), where Ny is the same Gaussian assumed in 
equation (25) and Ng is a Gaussian distribution (Lambas, Maddox, & Loveday 1992; Huizinga & 
van Albada 1992; Fasano et al. 1993; Ryden 2004): 


[( 6 -^- 

[ 2u2 J 


de . 


N^{e)de oc exp 


(26) 
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The function N{'y,€) is nonzero inside the region 0 < 7 < 1 and 0<e<l — 7 . A complication 
of fitting nonaxisymmetric models is the increased difficulty of deriving an inclination-corrected 
subset of spiral galaxies. If a galaxy is not an axisymmetric disk, it is impossible from photometry 
alone to determine absolutely whether or not its face-on semimajor axis R^q is greater than the 
cutoff value of 60". To illustrate the effects of nonaxisymmetry, consider a mildly triaxial disk, 
with e <C 1. For such a disk, P = l — e, D^l + 2ecos^ 4>, and the face-on semimajor axis is, from 
equation (17), 


R 


O 

20 


(1 -I- ecos^ 4>)R20 


1 -F {a/R2o)^^"{lnq 


e cos 2(j)) 



(27) 


Thus, the face-on isophotal radius R^q of a triaxial disk may be either smaller or larger than that 
of an axisymmetric disk with the same observed values of R 20 , ot, n, and q. Since R^q can be 
greater than the observed R 20 for a nonaxisymmetric disk, a sample complete to R 20 = 60'' is not 
guaranteed to be complete to R^q = 60". However, for mildly eccentric disks (e <C 1), the potential 
incompleteness effects should be small. 


To make the inclination correction, I begin by assuming values for the parameters /Zg, 

and (Te. I randomly select a value of 7 and of e from the appropriate distribution functions, and 
randomly choose a viewing angle {0,<p). I then compute (5 = 1 — e, using equation (7), 

E{9,(l)), using equation ( 8 ), and the apparent axis ratio q = \JDjE. I then have one possible 
(/9, D) pair corresponding to the observable apparent axis ratio q. Having binned the entire range 
of g, from g = 0 to g = 1, into 100 equal width bins, I can then place the {P,D) pair into the 
appropriate bin. I repeat the process of randomly choosing 7, e, 6, and (j), computing /3, D, E, and 
g, and putting the resulting (/3, D) pair into the appropriate bin in q, until I have done a total of 
~ 5 X 10® iterations of this process, keeping a maximum of 1000 (/ 3 , D) pairs for each bin in q. Now 
consider a specific galaxy with measured values of i? 20 ) ot, n, and qx- I go to the appropriate bin 
in axis ratio q and pull out {f3,D) pairs, determining how many of the pairs result in R^q > 60", 
using equation (17). A fraction E of the {P,D) pairs, where 0 < T < 1, will result in the face-on 
isophotal radius R^q of this galaxy being greater than the lower limit of 60". I then created binned 
distributions of qx and qs, with Ubin = 12, for comparison with model distributions. In the binned 
distributions of observed axis ratios, the axis ratio for each galaxy is given a weight equal to the 
fraction E. As a consequence, the total number n = ^ T of galaxies in the binned distribution 
may be non-integral. For convenience, I round the number n to the nearest integer by multiplying 
the number of galaxies per bin by the factor C = int(n -|- 0.5)/n. 

Finding the global best fit in the four-dimensional parameter space involves a 

tedious search. The search can be simplified by noting that for relatively thin, mildly elliptical 
disks, the values of and are determined by looking at galaxies that are close to edge-on, 
with q < 0.5; by contrast, the values of pe and are fixed by looking at galaxies that are close to 
face-on, with q > 0.5. Thus, the best-fitting values of and that I found by assuming e = 0 
should not be greatly changed when 0 < e <C 1. I therefore fix the values of and a-y for each 
sample of galaxies to be equal to those determined in the previous subsection, where galaxies were 
assumed to be axisymmetric. This reduces the parameter space to two dimensions. After fixing 
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and (T-y, I choose values of and iTc. I then randomly select n values of 7 from the distribution of 
equation (25), n values of e from the distribution of equation (27), and n random viewing angles. 
The apparent axis ratio (7(7, e, 9, (p) for each of n mock galaxies is computed from the formulas given 
in §2. The model axis ratios are then binned into Ubin = 12 equal-width bins. After repeating this 
procedure ~ 10^ times for a given (^i^, cjg) pair, I calculate the mean and standard deviation of the 
expected number of mock galaxies in each bin, and compute a probability for that particular 
set of parameters. 

Figure 7 shows the isoprobability contours in parameter space. The fit to the Kg band 

data is shown by the red lines; the fit to the B band data is shown by the blue lines. In each 
band, the best fitting values of and Ue in the displayed region of parameter space are indicated 
by an asterisk. The best fitting and fie for each galaxy sample, at each wavelength, are given 
in Table 1. The upper panel of Figure 7 displays the results for the complete sample of spirals, 
weighted by the appropriate fraction F for each galaxy. The best fits have « 0, and = 0.16 
[B band) and = 0.26 {Kg band). The most probable ellipticity for these skewed distributions is 
e ~ 0, but there is a long tail to higher ellipticities. The results for the 2MASS LGA spirals are 
in good agreement with those of Lambas, Maddox, & Loveday (1992), Fasano et al. (1993), and 
Ryden (2004), who find good fits for Gaussian distributions with = 0, and ~ 0.2. (I only 
computed solutions with //g > 0. It is mathematically permissible to have //e < 0. In the limit 
CTg, this leads to an exponential distribution function: oc exp(—e//ie), with a scale length 

/le = a‘l/\^f\ ^ fe- Exponential functions peaking at e = 0 do not give a significantly better fit 
than Gaussian functions peaking at e = 0.) 

The middle panel of Figure 7 shows the results of fitting a Gaussian distribution of ellipticities 
to the late spirals (Sc and later). This plot emphasizes that the apparent shapes of late spirals are 
consistent with axisymmetry. The best fit in each band has //g < 0.07 and = 0.02, but = 0, 
cJe = 0 also gives an excellent fit, with P > 0.5. Although perfect axisymmetry is not compulsory, 
there are constraints on how large the disk ellipticities can be. The region in parameter space for 
which P > 0.5 (demarcated by the heavy solid lines in the middle panel of Figure 7) is confined 
within the limits //g < 0.10, < 0.15. The parameters found by Fasano et al. (1993) for their 

sample of late-type galaxies, //e = 0 and = 0.09, are in good agreement with our results. However, 
the parameters found by Huizinga &: van Albada (1992) for their sample of Sc galaxies, //g = 0.12 
and (Te = 0.11, imply a greater mean ellipticity than is found in the 2MASS LGA sample. 

The bottom panel of Figure 7 shows the results of fitting a Gaussian distribution of ellipticities 
to the early spirals (Sbc and earlier). The best-fitting values of indicate a broad range of 
thicknesses: = 0.22 in the B band, and = 0.47 in the Kg band. Although late spirals 

have excellent fits with = 0, getting a good fit to the observed shapes of early spirals requires 
(Te > 0. With a broad tail stretching toward e = 1 — 7 (indicating a prolate shape), the best 
fit for early spirals in the Kg band doesn’t really represent a population of disks at all. If I 
quantify the triaxiality by the parameter T = (1 — /9^)/(l — 7^), the best fit in the Kg band, 
= (0.29,0.08,0.02,0.27), produces a mean triaxiality parameter (T) = 0.51, closer 
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to a prolate shape (T = 1) than an oblate shape (T = 0). The best fit to early spirals in the B 
band also has a significant mean triaxiality, (T) = 0.37. For comparison, the mean triaxiality for 
the best fit to late spirals is (T) = 0.04 in the Kg band and (T) = 0.13 in the B band. 


4.2. Lognormal parametric fits 

Satisfactory fits to the distribution of apparent axis ratios can also be found using a lognormal 
distribution of ellipticities, rather than a Gaussian. If I define rj = In e, then a lognormal distribution 
can be written in the form 

N^{e)de oc exp ^ dt] 

— (-^)f 

Motivation for the lognormal fit is provided by Andersen et al. (2001) and Andersen Sz Bershady 
(2002), who find, by combining photometric and kinematic information for nearly face-on spirals, 
that a lognormal distribution provides a good fit to the observed disk ellipticities. Finding the best 
values of firj and for each sample of spirals is done using the technique outlined in the previous 
subsection. Once again, the values of and used were those that gave the best ht in the 
axisymmetric case. 

Figure 8 displays the isoprobability contours in (/r^,cr^) space. For all subsamples and both 
colors, the lognormal distribution of ellipticities (equation 28) can provide as good a fit as the 
Gaussian distribution (equation 26). The best-fitting parameters and arj, as measured by a 
fit to the binned data, are indicated by asterisks in Figure 8 and are listed in Table 1. The best 
lognormal fit to each subsample yields a median ellipticity that is within ~ 15% of that provided by 
the best Gaussian fit. For all spirals, the median ellipticity is e ~ 0.10 in the B band and e ~ 0.16 
in the Kg band. 

The top panel of Figure 8 shows the results for the complete sample of spirals. The best fit in 
the B band, = —2.36 and arj = 0.96, is seen to be similar to the best fit found by Ryden (2004) 
for a sample of exponential galaxies in the Sloan Digital Sky Survey {firj = —2.06 and arj = 0.83 in 
the i band). It is also in agreement with the best fit to the sample of Andersen & Bershady (2002), 
corrected for their selection criteria (/r^ = —2.29 and = 1.04). The best fit in the Kg band, 
Hrj = —0.14 and arj = 2.74, requires a much broader spread in ellipticities in addition to a greater 
median ellipticity. However, the values of and arj are not very strongly constrained; excellent 
fits, with P > 0.5, are produced over a wide range of parameter space with firj ^ and u,; > 1. 

The middle panel of Figure 8 shows the results for the late spirals. The best-htting lognormal 
distribution in each band implies a small amount of disk ellipticity; the median ellipticity of late 
spirals is e ~ 0.07 in the B band and e ~ 0.02 in the K band. However, the shapes of late spirals 
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are seen, once again, to be consistent with axisymmetry —> oo, <7,^ = 0). The only region of 

parameter space that is strongly excluded is the lower right corner of the panel, where is large 
and cTr^ is small. Finally, the bottom panel of Figure 8 shows the results for the early spirals. The 
best fit in the B band implies a median ellipticity of e ~ 0.18 for the early spirals. The best fit in 
the Kg band implies a still greater degree of disk ellipticity, with a median ellipticity of e ~ 0.30. 


5. Implications and Discussion 

The 2MASS Large Galaxy Atlas provides a sample of nearby spiral galaxies with large angular 
size (i? 3 cr > 60'' at A ~ 2 microns). In this paper, I have used this relatively small sample of very 
well-resolved galaxies to estimate the distribution of apparent and intrinsic shapes of spiral galaxies. 
An intriguing result of this study is the difference in disk ellipticity e between disk galaxies of early 
and late Hubble type. 

If the entire population of 2MASS LGA spirals is considered, summed over all Hubble types, 
the distribution of ellipticities in the B band is consistent with earlier photometric studies using 
larger sample sizes (Lambas, Maddox, & Loveday 1992; Fasano et al. 1993; Ryden 2004). The Kg 
band ellipticity is greater, but the shape of the 2>a isophote in the Kg band can be influenced by 
light from triaxial bulges. The distribution N^[e) for shapes in the B band can be modeled as a 
Gaussian peaking at Ri 0 and width cje 0.16 or as a lognormal distribution with ~ —2.4 
and (7^ Ri 1.0; both these distributions imply a mean ellipticity e 0.1. When the late spirals, 
consisting mainly of Hubble type Sc, are examined, they are found to be perfectly consistent with 
axisymmetry in both the B and the Kg band. The early spirals, consisting mainly of types Sb and 
Sbc, show a high degree of ellipticity, even in the B band, in which the mean ellipticity is e ~ 0.18. 

It has been noted by Franx &: de Zeeuw (1992) that a disk ellipticity of e > 0.1 cannot be 
reconciled with the observed small scatter in the Tully-Fisher relation (Tully & Fisher 1977) if 
the disk ellipticity is assumed to trace the potential ellipticity. Suppose that stars and gas are 
on closed orbits in a logarithmic potential with rotation speed Vc- If the potential has a small 
ellipticity <C 1 in the orbital plane, then the integrated line profile from the disk will have a 
width W = 2uc(l — cos 2(/>) sin 0 when viewed from a position angle (0,</>) (Franx & de Zeeuw 
1992). The difference in line width from that produced in a purely circular disk will produce a 
scatter in the observed Tully-Fisher relation. If all disks have = 0.1, the expected scatter is 0.3 
mag, even if the inclination has been determined accurately from kinematic information (Franx 
&: de Zeeuw 1992). In most studies of the Tully-Fisher relation, the inclination is determined 
photometrically, from the apparent axis ratio of the disk, assuming (perhaps erroneously) that the 
disk is axisymmetric. If inclinations are determined in this way, the scatter in the Tully-Fisher 
relation will be even greater. 

If galaxies had infinitesimally thin circular disks, the relation between apparent axis ratio q 
and inclination 6 would be = cos 9. If the thin disks actually have ellipticity e and are viewed 
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at high inclination (sin^ 0 ^ 2e), the apparent axis ratio will be g ~ cos 0(1 — ecos2(^). Thus, 
if the inclination is estimated by the relation 0 phot = cos“^ g, a fractional error oc ecos 2 i;^ will 
be introduced into cos0phot and an error oc ecot^0cos2(^ will be introduced into sin0phot- Since 
the observed velocity width W must be divided by sin0 to find the rotation speed Vc, the error 
produced by using sin0phot instead of sin0 will be negligible only when cot^ 0 <C 1. Thus, most 
studies of the Tully-Fisher relation use galaxies of high inclination (0phot ^ 45°). The Tully-Fisher 
relation for nearly face-on galaxies (see, for instance, Andersen & Bershady (2003)) can only be 
determined if the inclinations are determined kinematically. 

In Figures 7 and 8 , the dashed lines indicate the expected amount of scatter in the Tully- 
Fisher relation if the potential ellipticity is given by either a Gaussian distribution (Figure 7) or by 
a lognormal distribution (Figure 8 ). Starting at the lower left corner of each panel, the contours 
are drawn at the levels 0.25 mag, 0.5 mag, 0.75 mag, and 1.0 mag. The inclinations are assumed 
to be estimated from the apparent axis ratio. Since the disks are not always infinitesimally thin, 
the inclination 0 phot is estimated using the relation 

cos^ 0phot = 

J- 

Ultrathin galaxies with q < are assumed to have cos0 = 1. All galaxies with 0phot < 45° are 
discarded. 

If all spirals are considered together, the best fitting ellipticity distributions - either Gaussian 
or lognormal - imply approximately 1 mag of scatter if the potential ellipticity equals the disk 
ellipticity: this can be seen in the upper panels of Figures 7 and 8 . The best fitting ellipticity for 
the late spirals, as seen in the middle panels of Figures 7 and 8 , implies only 0.3 mag of scatter. 
Finally, the best fitting ellipticity for the early spirals, seen in the lower panels of Figures 7 and 
8 , leads to 1.4 mag of scatter in the Tully-Fisher relation. By comparing the predicted scatter 
from disk ellipticity to the actual scatter in the Tully-Fisher relation for these galaxies, I can place 
constraints on how much of the scatter can be attributed to the disk ellipticity. 

To create a Tully-Fisher plot for the spirals in the 2MASS LGA sample, I took the relevant 
astrophysical parameters from the HyperLeda database.^ As a measure of galaxy luminosity, I 
used the I band absolute magnitude; multiband studies of the Tully-Fisher relation indicate that 
the observed scatter is minimized in or near the near-infrared. The I band absolute magnitude 
is computed from the apparent magnitude mj (corrected for galactic extinction and internal ex¬ 
tinction) and the distance modulus m — M. The distance modulus was computed from the radial 
velocity uyir corrected for infall of the Local Group toward Virgo, assuming a Hubble constant 
Hq = 70 km s“^ Mpc~^; for galaxies with uvir < 500kms“^, the distance modulus was taken from 
the literature, with preference given to Cepheid distances. As a measure of rotation speed, I 
used the maximum rotation velocity Vm determined from the 21 cm line of neutral hydrogen. In 


^See http://leda.univ-lyonl.fr/ 
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computing Vm, the inclination used is the photometric estimate 0phot) using the axis ratio of the 
25 mag arcsec“^ isophote as the apparent axis ratio q. 

Figure 9 shows Mj as a function of log Vm for the 2MASS LGA spirals with > 60", assuming 
axisymmetry; galaxies with 0phot < 45° are omitted, as are galaxies without I band photometry 
in HyperLeda. A total of n = 128 galaxies meet all these criteria. The solid line in Figure 9 is 
the best fitting Tully-Fisher relation, found using an unweighted inverse fit (Schechter 1980; Tully 
Sz Pierce 2000; Kannappan et al. 2002). The best Tully-Fisher relation for the sample of n = 128 
galaxies is 

M,(all) = (-22.43 ± 0.09) - (9.42 ± 0.77) log^o > (30) 

with a scatter in Mj of crrms = 0.94 mag. If only the n = 37 late spirals are considered, an 
unweighted inverse fit yields 

M,(late) = (-22.54 ± 0.15) - (9.42 ± 0.85) logip > (31) 

with a scatter in Mj of Urms = 0.62 mag. Finally, if only the n = 91 early spirals are considered, 
the fit is 

A7,(early) = (-22.40 ± 0.12) - (9.97 ± 1.38) logip > (32) 

with a scatter in Mj of Urms = 1-10 mag. The late and early type spirals have Tully-Fisher relations 
with statistically indistinguishable slopes; however, the early spirals have a much greater scatter. 

The observed scatter, Urms, is not due entirely to the nonaxisymmetry of disks. In addition to 
an intrinsic scatter fjeii, due to the ellipticity of the disk and the potential, there is also a contribution 
(Terr due to errors in observation and interpretation. For instance, the absolute magnitude M/ is 
subject to errors in flux measurement, errors in determining distance, and errors in extinction 
correction. The rotation speed Vm is subject to errors in measuring 21 cm line widths and errors in 
converting line widths to deduced rotation speeds. Let us suppose that the two sources of scatter 
add in quadrature: = (T^^ -|- If (Terr is the same for both early and late spirals, then 

(Teu(early)^ - fTeu(late)^ = fTrms(early)^ - fTrms(late)^ = (0.91 mag)^ . (33) 

If late spirals are perfectly axisymmetric, which is permitted by the data, then (Teii(late) = Omag 
and (Tell(early) = 0-91 mag. The largest possible values of (Teu come if I naively assume that (Terr = 0; 
in this limiting case, (Ten(late) = 0.62 mag and (Teu(early) = 1.10 mag. 

For late spirals, the Tully-Fisher relation thus implies that the intrinsic scatter due to ellipticity 
lies in the range Ueii = 0 —> 0.62 mag. This is consistent with the value (Teii = 0.3 mag derived from 
the best fitting distributions of disk ellipticity. For early spirals, the Tully-Fisher relation implies 
an intrinsic scatter due to ellipticity in the range (Teu = 0.91 —> 1.10 mag. This is smaller than 
the value (Teii = 1.4 mag derived from the best fitting distributions of disk ellipticity. However, in 
the lower panels of Figures 7 and 8, the band where (Tell = 0.91 —> 1.10 overlaps the region where 
P > 0.1; thus, the ellipticities derived from the apparent shapes and from the Tully-Fisher relation 
are not discrepant at a high confidence level. 
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The difference in disk ellipticity between earlier and later spirals, and the resulting difference 
in Tully-Fisher scatter, helps to explain the apparent discrepancy between the relatively large disk 
ellipticity (e ~ 0.1) seen in photometric studies (Binney & de Vaucouleurs 1981; Grosbpl 1985; 
Lambas, Maddox, & Loveday 1992; Fasano et al. 1993; Alam & Ryden 2002; Ryden 2004) and the 
relatively small intrinsic scatter (cr < 0.3 mag) seen in some studies of the Tully-Fisher relation. 
For instance, Verheijen (2001), in his study of spirals in the Ursa Major cluster, found that his DE 
(distance estimator) subsample of spirals had a Tully-Fisher relation with 0.26 magnitudes of scatter 
in the I band; the intrinsic scatter had a most likely value of only 0.06 magnitudes. This small 
scatter was partly due to the strict selection criteria imposed; galaxies were excluded if they were 
obviously interacting, had prominent bars, or had irregular outer isophotes. Galaxies of Hubble 
types earlier than Sb or later than Sd were also excluded. The scatter was further minimized by 
the fact that the resulting subsample of 16 spirals contained 10 galaxies of type Sc, Scd, and Sd 
(late spirals, in my classification) and only 6 galaxies of type Sb and Sbc. If the spirals in the 
2MASS LGA are characteristic of all spiral galaxies, then scatter in the Tully-Fisher relation can 
be reduced by excluding spirals of Hubble type earlier than Sc. 

In addition to the practical implications for minimizing scatter in the Tully-Fisher relation, the 
difference in ellipticity between early and late spirals provides an intriguing clue for the formation 
and evolution of spiral galaxies. Why should the Sb and Sbc galaxies dominating the early spiral 
sample be more elliptical at R 25 than the Sc galaxies? One possibility is that the triaxiality of 
a dark halo’s potential affects the Hubble type of the spiral embedded within it. Kormendy & 
Kennicutt (2004) estimate that the majority of spirals with Hubble type of Sb and later contain 
“pseudobulges” rather than classical bulges. Pseudobulges are formed as gas is transported to 
small radii by nonaxisymmetric structures such as triaxial halos, spiral structure, elliptical disks, 
and bars. Simulations of disk galaxies in triaxial dark halos (El-Zant & Hassler 1998; El-Zant 2001) 
indicate that pseudobulges grow by secular evolution on a timescale of a few Gyr. The Hubble type 
of the resulting galaxy depends on both the halo core radius and the potential asymmetry, with 
greater halo nonaxisymmetry leading to larger bulges. 

High resolution n-body simulations of virialized dark halos in the mass range 10^^ Mq < M < 
10^“^ Mq indicate that virialized halos are well described as triaxial shapes (Jing Sz Suto 2002; 
Allgood et al. 2005). More massive halos are further from spherical, on average. In a lambda GDM 
simulation, the mean short-to-long axis ratio of a halo with mass M at the present day {z = 0) 
is ( 7 ) Ri , where M* = 1.2 x 10^^ Mq is the characteristic nonlinear mass scale 

(Allgood et al. 2005). The mean intermediate-to-long axis ratio also decreases with increasing mass, 
going from (/3) Ri 0.76 at M = O.IM* to (/3) Ri 0.71 at M = M*. (Dissipation by baryons tends 
to compress the dark halo into a more nearly spherical shape in its central regions (Kazantzidis 
et al. 2004), but will not destroy the trend that lower-mass halos are rounder on average in the 
plane of the baryonic disk.) At a given halo mass, there exists a fairly wide spread in f5 around 
the mean value, with standard deviation ~ 0.1 for (/3) ~ 0.7 (Allgood et al. 2005). The nearly 
axisymmetric halos will give rise to nearly circular disks and smaller pseudobulges, such as those 
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seen in Sc galaxies. The more nonaxisymmetric halos will give rise to elliptical disks and larger 
pseudobulges, such as those seen in Sb and Sbc galaxies. The mean mass of Sb spirals is greater 
than that of Sc spirals; however, the spread in masses for a given Hubble type is greater than the 
difference in the mean (Roberts Sz Haynes 1994). This would be expected if the difference between 
early (Sb and Sbc) spirals and late (Sc) spirals is not due primarily to differences in halo mass, but 
rather is due to differences in halo nonaxisymmetry. 


This publication makes use of data products from the Two Micron All Sky Survey, which is a 
joint project of the University of Massachusetts and the Infrared Processing and Analysis Center / 
California Institute of Technology, funded by the National Aeronautics and Space Administration 
and the National Science Foundation. It also makes use of the HyperLeda galaxy database. James 
Pizagno provided useful Tully-Fisher insights. 
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Table 1. Best Parametric Fits 


Sample 

Band 


CT/y 


CTe 


CTj} 

All 

Ks 

0.26 

0.09 

0.01 

0.26 

-0.14 

2.74 


B 

0.17 

0.03 

0.01 

0.16 

-2.36 

0.96 

Late 

Ks 

0.18 

0.06 

0.02 

0.02 

-3.86 

0.74 


B 

0.14 

0.02 

0.07 

0.02 

-2.63 

0.13 

Early 

Ks 

0.29 

0.08 

0.02 

0.47 

-0.14 

1.74 


B 

0.20 

0.04 

0.09 

0.22 

-1.77 

0.83 
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Fig. 1.— Axis ratio of the 25 mag/arcsec^ isophote in the B band versus the axis ratio of the 3cr 
isophote in the Kg band. Squares indicate late spirals (Sc or later); triangles indicate early spirals. 
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Fig. 2.— The semimajor axis of the 3a isophote, measured in units of the half-light radius Re 
of the best Sersic profile, plotted as a function of the axis ratio qk- Squares indicate late spirals; 
the dashed line is the best fit to their R^f^/Re vs. In(gi^) relation. Triangles indicate early spirals; 
the dotted line is the best fit to their R^er/Re vs. ln{qK) relation. The labeled outlier, NGC 1068, 
is excluded from the fit. 
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Fig. 3.— Upper: Distribution of apparent axis ratios for the inclination-adjusted sample of 2MASS 
LGA spirals. Red lines indicate the Kg axis ratios; blue lines indicate the B axis ratios. The 
solid line in each case is the best ht, while the dotted lines indicate the 80% confidence interval 
estimated from bootstrap resampling. The horizontal error bars show the initial kernel width h. 
Lower: Distribution of intrinsic axis ratios, assuming the spirals are randomly oriented oblate 
spheroids. 
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Fig. 4.— The same as Figure 3, including only the n = 63 spirals of Hubble type Sc and later. 
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Fig. 5.— The same as Figure 3, including only the n = 130 spirals of Hubble type Sbc and earlier. 
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Fig. 6.— Top: Asterisks indicate the best fitting values of and to the complete sample of 
LGA spirals; the heavy solid line, light solid line, and dotted line indicate the P = 0.5, P = 0.1, 
and P = 0.01 isoprobability contours, found by a fit- Red = shapes deduced from Kg data; blue 
= shapes deduced from B data. Middle: Same, but for spirals of Hubble type Sc or later. Bottom: 
Same, but for spirals of type Sbc or earlier. ‘F’ = fits of Fasano et al. (1993). ‘H’ = fit of Huizinga 
&: van Albada (1992) for Sc spirals. 
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Fig. 7.— Top: Asterisks indicate the best fitting values of /i^ and to the complete sample of 
LGA spirals. The heavy solid line, light solid line, and dotted line indicate the P = 0.5, P = 0.1, 
and P = 0.01 isoprobability contours. Red = shapes from Kg data; blue = shapes from B data. 
Dashed lines show scatter in the Tully-Fisher relation (see §5); contours are drawn, starting at the 
lower left, at 0.25, 0.5, 0.75, and 1.0 magnitude. Middle: Same, but for spirals of Hubble type Sc 
or later. Bottom: Same, but for spirals of type Sbc or earlier. ‘F’ = Fasano et al. (1993); ‘L’ = 
Lambas, Maddox, & Loveday (1992); ‘R’ = Ryden (2004); ‘H’ = Huizinga & van Albada (1992). 
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Fig. 8.— Top: Asterisks indicate the best fitting values of p, and cr to the complete sample of 
LGA spirals, assuming a lognormal distribution of disk ellipticities. The heavy solid line, light 
solid line, and dotted line indicate the P = 0.5, P = 0.1, and P = 0.01 isoprobability contours. 
Red = shapes deduced from Kg data; blue = shapes from B data. Dashed lines show the scatter 
in the Tully-Fisher relation; contours are drawn, starting at the lower left, at 0.25, 0.5, 0.75, and 
1.0 magnitude. Middle: Same, but for spirals of Hubble type Sc or later. Bottom: Same, but for 
spirals of type Sbc or earlier. ‘A’ = best fit to data of Andersen et al. (2001); ‘R’ = best fit of 
Ryden (2004). 
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Fig. 9.— The corrected I band absolute magnitude versus the logarithm of the maximum 21 cm 
rotation velocity. Squares indicate late spirals; triangles indicate early spirals. The solid line is the 
inverse unweighted fit to all the points (n = 128), the dashed line is the fit to the late spirals only 
(n = 37), and the dotted line is the fit to the early spirals (n = 91). 





